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INTERACTIONS OF THE MAGNETOSPHERES OF STARS AND CLOSE-IN GIANT PLANETS 
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ABSTRACT 

Since the first discovery of an extrasolar planetary system more than a decade ago, hundreds more 
have been discovered. Surprisingly, many of these systems harbor Jupiter-class gas giants located 
close to the central star, at distances of 0.1 AU or less. Observations of chromospheric 'hot spots' 
that rotate in phase with the planetary orbit, and elevated stellar X-ray luminosities, suggest that 
these close-in planets significantly affect the structure of the outer atmosphere of the star through 
interactions between the stellar magnetic field and the planetary magnetosphere. Here we carry out 
the first detailed three-dimensional MagnetoHydroHynamics (MHD) simulation containing the two 
magnetic bodies and explore the consequences of such interactions on the steady-state coronal struc- 
ture. The simulations reproduce the observable features of 1) increase in the total X-ray luminosity, 2) 
appearance of coronal hot spots, and 3) phase shift of these spots with respect to the direction of the 
planet. The proximate cause of these is an increase in the density of coronal plasma in the direction 
of the planet, which prevents the corona from expanding and leaking away this plasma via a stellar 
wind. The simulations produce significant low temperature heating. By including dynamical effects, 
such as the planetary orbital motion, the simulation should better reproduce the observed coronal 
heating. 

Subject headings: Stars: coronae — (Stars:) Planetary Systems — Stars: Activity 



1. INTRODUCTION 

The structure and heating of the solar corona, as well 
as the acceleration of the solar wind, are influenced 
by the structure and topology of the large-scale coro- 
nal magnetic field. On this basis, t he existence of a 
plane t at a distance of .1 AU or less (|Mavor fc Quelod 
119951 : iMayor et al.l l2003h . with a strong internal mag- 
netic field is expected to have a significant effect on the 
stellar magnetosphe re, which is contr olled by the mag- 
netic field structure (jCuntz et al.ll2000i ). In recent years, 
some signatures of this star-planet i nteraction (SPI) have 
been observed. Shkolnik et al. ("Shkol nik et~all 120031 . 
[^05a,b, 2008(1 have reported on modulations in the Ca II 
K emission line, an indicator for chromospheric activity. 
They find enhancements in the line intensity that have 
the same period as the planetary orbital motion, though 
sometimes with a significant non-zero phase-shift. The 
cause is deemed magnetic and not tidal because of the 
lack of an equivalent hot spot offset in phase by 180°. 
In addition, a statistical survey of the X-ray fiuxes from 
stars with close-in planets has found them enhanced by 
a 30-400% on average over typical flux es from similar 
stars with planets that are not close-in (|Kashvap et al.l 
I2008D Direct X-rav o bservations of the HD 179949 sys- 
tem (jSaar et al.ll2008f l showed that the SPI contributed 
«30% to the emission at a mean temperature of « 1 keV. 

Some analytical and semi-empirical arguments have 
been advanced to explain these observations. One posits 
that particles are accelerated along magnetic field lines 
that connect the star and planet, creati ng hot spots 
where they hit the chromospheric layer ( Cuntz et al.l 
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120001: iLanzal 120081: iCranmer fc SaaH 120071) . As a result, 
hot spots are observed generally in phase with the plan- 
etary orbit, but with the capacity to have large offsets, 
depending on the exact structure of the magnetic field 
between the star and planet. Another shows that tran- 
sition of field lines from a high-helicity state to a linear 
force-free state is e nergetically a dequate to power the en- 
hanced intensities (|Lanzal[2009f) . The detailed behavior 
of the dynamical interaction of coronal and wind plasma 
with two magnetic field systems is, however, very difficult 
to realize with idealized models. The problem properly 
requires simultaneous descriptions of both the stellar and 
the planetary magnetospheres, the planetary orbital mo- 
tion, and often asynchronous stellar rotation, together 
with a self-consistent stellar wind solution. 

Here we describe an initial simulation of the magnetic 
star-planet interaction. We use idealized test cases to 
study the fundamental changes in the steady-state coro- 
nal structure due to the presence of the planet and its 
magnetic field. The dynamical interaction due to the 
planetary orbital motion is captured in an indirect man- 



2. SIMULATION 

The numerical simulation has been performed using 
the University of M ichigan Solar Corona (SC) model 
(jCohen et al.l I2007D . w hich is based on t he BATS-R- 
US global MHD code (jPowell et al.l 119990 and is part 
of the Space We ather Modeling Framework (SWMF) 
(|T6th et al.ll2005l ). The model solves the set of magneto- 
hydrodynamic equations on a Cartesian grid using adap- 
tive mesh refinement (AMR) technology. This model 
has been extensively validated for the solar corona us- 
ing coronal observations and in-situ solar wind measure- 
ments taken at 1 AU (,Cohen et al.|[200a) . We assume 
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that the particular physical description of the coronal 
heating and wind acceleration is not crucial to study the 
change in the existing coronal structure due to the planet. 
It is important to mention that we use a global model for 
the corona that cannot reproduce realistic chromospheric 
emission due to heating of coronal loops. We also do 
not fully describe the observed coronal heating, since for 
example no input from magnetic reconnection or loop 
footpoint motion is included. Thus, while we adopt the 
physical parameters of some real systems in the model- 
ing, we do not expect the models to fully reproduce all 
aspects of observations (in particular, details of the tem- 
perature and level of the emissions) at this point. The 
full physical des cription of the model and its limitations 
can be found in ICohen et all ([2007I . [200l 

We performed several different numerical simulations, 
of which we highlight two here. 'Case A': both the stel- 
lar and planetary magnetic fields are perfectly aligned 
dipoles. We set the stellar polar field to be 5G and 
the planetary polar field to be antiparallel at —2G (i.e., 
opposite to the stellar dipole). The planetary mag- 
netic field is weaker than Jupiter's, and follows the as- 
sumption that hot-Jupiters are expected (but not re- 
quired) to have lower spin rates d ue to tidal locking, and 
thus have weaker magnetic fields (lSanchez-Lavega|l2005 
iGriefimeier et al.l 12004 lOlson fc ChristensenI 120061 ). We 
note that a simulation in which the planetary dipole was 
set to be in the same direction with the stellar dipole 
resulted in a quantitatively similar solution as in this 
case. 'Case B': the planetary magnetic field is a perfect 
dipole and the stellar magnetic field is driven by solar 
magnetic synoptic map (magnetogram). This map con- 
tains measurements of the photospheric radial magnetic 
field taken during solar maximum (Carrington Rotation 
CR2010, very active Sun). The use of a magnetic syn- 
optic map enables us to generate a realistic, Sun-like, 
three-dimensional magnetic field. 

In Case A, we mimic the relative motion between the 
planet and the background plasma by fixing the planet 
and rotating the star and the coronal plasma in the in- 
ertial frame. This way, the planet orbits the star back- 
wards in the frame rotating with the star. This is done 
due to the fact that the actual orbital motion of the 
planet requires time-dependent boundary conditions. We 
plan to implement this technical improvement in future 
simulations. For the sake of definiteness, we partially 
match the parameters o f the system to the observed pa- 
rameters of HD 179949 (|Mavor et al. 2003), which is an 
F8V type star. We use the following stellar parameters: 
Mi, = 1.28 Mq, i?^ = 1.19 i?0, and stellar rotation pe- 
riod of 3 d. In the HD 179949 system, the planet is lo- 
cated at a distance of 0.045 AU (9.65 stellar radii), and a 
60° phase-lead of a chromospheric hot spot is observed. 
The planetary parameter Msini = 0.98 Mj has not been 
used here. In Case B, we fix the planet relative to the 
star and run the simulation in the frame rotating with 
the star. In this case we use solar parameters except for 
the planetary properties, which are the same as in Case 
A. This case represents a steady-state, large scale inter- 
action of a Sun-like star w ith with a tidally - locked extra- 
solar planet (e.g., r Boo; (|Saar et al.ll2004l : iCatala et al.l 
I2007D ) located at the same distance as before. 

In both simulations, the boundary condition for the 
planetary plasma number density and temperature were 



no = 10^° cm and Tq = lO"' K respectively 
(|Murrav-Clav et ani2009D . The stellar boundary condi- 
tions were no = 10^ cm~^ and Tq = 3.5 • 10^ K respec- 
tively, based on pre vious simulation of the solar corona 
(jCohen et al.l 120071 ) . To further aide interpretation of 
the results, we performed two additional simulations as 
a reference, identical to the cases above, but with the 
planet removed; i.e., considering the star with just the 
5G dipolar field, and with the CR2010 magnetogram. 

In each simulation, the set of MHD equations is solved 
until convergence. The end result is a three-dimensional, 
steady-state solution for the particular system that in- 
cludes all the MHD variables (density, pressure, velocity, 
and magnetic field). Since the MHD solution contains 
the values for n^ and T at each spatial cell, we can per- 
form the line-of-sight integration to obtain the predicted 
X-ray emissions for a particular view angle. The integra- 
tion takes in to account cells in front of the star but omits 
cells behind it. We repeat this procedure for different 
view angles to mimic the predicted X-ray flux as the sys- 
tem rotates. The X-ray flux, f^, is calculated as the line- 
of-sight integral fx ^ / '^e ' P{T)dl. Here dl is the line- 
of-sight depth and we have used a piecewise linear ap- 
proximation to the radiative los s P(T) for a plasma wi th 
solar photospheric abundances (jCiampapa et al.lll996[ ). 

3. RESULTS 

Simulation results are illustrated in Figure [H Top and 
middle panels show the three-dimensional solutions ex- 
cluding and including the planet, respectively. Bottom 
panels show the difference in X-ray flux for tempera- 
ture range of log{T) = 6.04 — 6.23 with and without the 
planet. Left panels show results of Case A, while right 
panels show results of Case B. 

Considering first Case A, we note a key difference be- 
tween simulations with and without the planet. In the 
former, magnetic field lines are conspicuously brought in 
toward the planet and are constricted due to the pres- 
ence of the planetary magnetosphere. This has a palpa- 
ble influence on the coronal electron number density, n^ , 
which now increases azimuthally approaching the star- 
planet line to form an X-ray 'bright spot' facing the 
planet. This solution is qualitatively very similar in Case 
B, where there is also a clear longitudinal concentration 
of the plasma density, and consequently the X-ray flux 
(which is proportional to Tig). While we cannot simu- 
late chromospheric emission with our current models, its 
surface intensity distribution on the Sun follows closely 
those regions of the disk that are brighter at EUV and 
X-ray wavelengths. The results of these simulations are, 
then, fully consistent with the observed location of chro- 
mospheric hot spots seen in-phase with the planetary 
orbit. 

The longitudinal brightening effect is also clear in Fig- 
ure [21 where the line-of-sight X-ray flux originating at 
different plasma temperatures is shown as a function of 
viewing angle for the simulations including planets. In 
the idealized dipole field case, the emission is much more 
intense when the corona is viewed from the direction of 
the planet than from the opposite direction when the 
star hides the brighter parts of the corona that form 
the hot spot in phase with the planet (a difference of 
25-35% in observed X-ray flux). The simulation based 
on the realistic, complex magnetic field results in a hot 
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spot (15-30% difference in X-ray flux) shifted by about 
60° relative to the star-planet Une. This suggests that 
the phase shifts between hot spots and planetary orbital 
phase seen in SPI observations are probably due to the 
complexity of the stellar coronal magnetic fields and the 
consequent complexity of the magnetic connectivity be- 
tween star and the planet. 

Also of interest is the magnitude of the X-ray flux en- 
hancement seen compared to the case where there is no 
planet. In the ideal dipole field case, we find enhance- 
ments of at least 10% in base emission and as much as 
«80% relative to the no-planet case, and a contrast be- 
tween minimum and maximum of 20-30%. These num- 
bers are consistent with the observed signatures of SPI 
found in X-ray observations. When solutions driven by 
realistic magnetogram data are considered, we find that 
enhancements in base emission by a factor of 10 are pos- 
sible. The simulation with realistic, complex magnetic 
field results in much higher density enhancement and 
closing of coronal loops compared to the dipolar case. 

The MHD model, however, does not take into account 
the detailed physics of smaller scale X-ray emitting coro- 
nal loops, and the magnitude of the predicted X-ray en- 
hancements should be considered approximate. The den- 
sity enhancement in the closed-field zone provides a vital 
medium which can be heated to produce the SPI effect. 
The compression alone produces significant increase in 
radiative loss at low temperatures; localized dynamical 
effects, not modeled here, such as megnetic reconnection 
due to the planetary orbital motion, and particle acceler- 
ation then likely further heat the plasma to the observed 
^ 1 keV temperatures (iSaar et al.ll2008[ ). 



4. CONCLUSIONS 

In conclusion, we find that a dominant physical effect 
creating observable time-variable spi signatures is that 
the existence of the planet and its magnetosphere, close 
to the star, prevents the expansion of the stellar coronal 
magnetic field and the acceleration of the stellar wind. 
The pressure gradient is not as large as it would be in the 
absence of the planet, so the coronal field lines that would 
be opened by the wind remain closed and the plasma in 
these loops do not escape. This effect alone reproduces 
three observable feature: 1) Enhancement of total X-ray 
flux, 2) Appearance of coronal hot spots, 3) Phase-shift 
of the hot spots from star-planet line. The density en- 
hancement results in low temperature coronal heating. 
We will further develop the model to include the plane- 
tary orbital motion in order to capture more dynamical 
effects of SPI. 
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Fig. 1. — Illustration of the MHD simulation results for Case A (left), and for Case B (right). The top panels show the reference 
simulations without the planet, while the simulations including the planet are shown in the middle panels (planetary surface is marked by 
black circle). Color contours represent plasma number density on the y=0 plane, and yellow three-dimensional streamlines represent stellar 
magnetic field lines that originate on the intersection of the y=0 plane with the stellar surface. The inner boundary of the simulation 
domain (the stellar surface) is represented by a red sphere, and the view is from an angle of 45° . In the middle panels, the magnetic field 
lines that connect to the planetary magnetosphere are drawn in pink. The bottom panels show the difference in X-ray emission for plasma 
in the logarithmic temperature range logT = 6.04 — 6.23 between the case with and without the planet. Black circles represent the stellar 
surface and indicate that the 'hot spots' are located in the low corona. 
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Fig. 2. — Integrated LOS X-ray flux for the dipolar case with relative orbital motion (top) and for the realistic solar maximum case 
(bottom) as a function of view angle. The planet is located at 90° and is hidden by the star at 270° . Left panels show the flux for three 
temperature bins, where each curve is normalized to its own maximum value. Right panels show the ratio of the X-ray flux with and 
without the planet for each temperature bin as a function of view angle. 



